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The availability of computational modeling tools for subatomic physics (Form, FeynArts, 
FormCalc, and FeynCalc) has made it possible to perform sophisticated calculations in per- 
turbative quantum field theory. We have adapted these packages in order to apply them 
to the effective chiral field theory of hadronic interactions. A detailed description of this 
Computational Hadronic Model is presented here, along with sample calculations. 

I. INTRODUCTION 

Chiral Perturbation Theory (ChPT) has been tremendously successful in describing low-energy 
hadronic interactions in the non-perturbative regime of QCD. It provides for a systematic and well-defined 
perturbative expansion incorporating interaction physics driven by the symmetries of QCD in kinematic 
regimes where scales are well-defined. 

Calculations have been performed for many problems to Next-to-Leading-Order (NLO) in a variety of 
bases, but the challenge is often in determining the range of valid and contributing degrees of freedom to 
a given problem. 

With recent developments in the automatization of the NLO calculations in perturbative field theory, 
it is reasonable to consider the possibility that these methods could also be applied to ChPT and thus 
allow for a broader application of the theory. 

To date, there are several packages available (FeynArts [H], FormCalc [3], FeynCalc [1] and Form jl]) 
allowing us to produce semi-automatic calculations in particle physics. These packages have been applied 
in various calculations relevant to studies of the Standard Model 

Although FeynArts and FormCalc were originally designed for Standard Model calculations, the flexi- 
bility of the programs allows us to extend them to interactions appropriate for the hadronic sector. This 
was a main reason to use FeynArts and FormCalc as a base languages for the automatization of chiral 
hadronic calculations. Hence the main purpose of this work is to develop an extension to FeynArts and 
FormCalc which will perform automatic analytical calculations in ChPT. 

We will begin by presenting a summary of the properties of Chiral Perturbation Theory, followed by 
a description of how the interactions can be incorporated into a computational hadronic model. We will 
then present some sample applications and tests of the model, along with plans for future development. 

II. EFFECTIVE CHIRAL LAGRANGIAN 

Let us start with a general review of the formalism we use to describe the dynamics of the strong 
interactions. Spontaneous symmetry breaking of SU (3)^ (8i SU due to the pseudo-Goldstone scalar 
bosons into SU (3)y can be described by the Lagrangian from 



= 

Eq.([T]) represents the first term of the effective Lagrangian from which is the only allowed term with 
two derivatives. And additional term Tr^^^ is a constant due to the equivalence theorem and is excluded 
from the effective Lagrangian. Here, ~ 135 MeV is the pion decay constant and ^ field is given by 

^ = g2iP/^^ (2) 
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where P is the pseudo-Goldstone boson octet: 



75^ + 



P=\ - 75^-75" ^° I- (3) 
Taking into account \J (1)q of electromagnetism, the covariant derivative is given by 

D^ = a^ + a;,[Q,...], (4) 

where A enters as the electromagnetic vector field potential and Q is a charge operator. Unlike the baryon 
field, chiral symmetry transformation for the pseudo-Goldstone boson field exists and is well defined as 
^ LJ2R^ ■ To introduce baryon field into the effective Lagrangian uniquely, the new field ^ is defined 
as = The chiral symmetry transformations for the field ^ can be determined in the a new basis U 
(see [HI) as 

C L^U^ = UiR\ (5) 

Here, the unitary matrix U is implicitly defined from Eq.(l5]) in terms of L, i2, and ^, and in this case the 
chiral transformation for the baryon field is unique, and chosen to be — > UBU^. Choice of the basis U 
is preferable for the effective Lagrangian with baryons, because in this basis pions are coupled through 
the derivative type coupling only. Consequently, the effective Lagrangian for the baryons can be written 
in terms of the vector fields (vector and axial- vector A^): 



A, = l (^D^^t _ ^tz)^^) = l.d,P + -Ip {d,P) P-^ {d,P, P'} 



(6) 



The leading order baryon Lagrangian is given by 

^bI = -iTrB^B + m^Tr BB + 2D Tr ^7^75 {A,,, B} + 2F Tr B-f^'-f5 [A^, B] , (7) 
where B is SU (3) octet of baryons given by 




with covariant derivative defined as ^ =d^ + [Km---]- The strong coupling constants {F, D} of the La- 
grangian from Eq.© have been determined in to be F = 0.40 ± 0.03 and D = 0.61 ± 0.04. At the 
lowest order of chiral perturbation theory, the strong interaction of the decuplet T^{J^ = 3/2=*=) can be 
described by the following Lagrangian: 

4'°^ = -if' pT^ + mrf^'T^ + C {BA^T^T^^ + f^A^T^B) + inf^Yl^^uT^T^,, (9) 
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J^=3/2+ J^=3/2~ 

where Fs = 1 for and = 75 for . The decuplet states are given by 



■^333 



ft' 



rpl33 _ J_^*0 2^233 _ 1 



2^113 



jilll _ ^++ j^n2 _ ^_^+ 2^122 _ rp222 _ 

The strong couphng constants \C\ = 1.2 and Ti ~ — 1.8 have been indirectly extracted from the experimental 



data by [13|| through the loop corrections of the strong decay of decuplet to octet of baryons in the 



framework of Heavy Baryon [HBxPT] . 

III. COMPUTATIONAL HADRONIC MODEL 

The effective chiral theory of the strong interactions formulated above can be applied towards a vast 
number of the hadronic reactions in the nonperturbative regime. However, in this formalism, to do 
calculations up to NLO by hand requires a tremendous amount of work. Because these calculations are 
rather algorithmic, it is only reasonable to raise the question whatever some degree of the automatization 
of the present calculations in the ChPT is currently possible. We address this question by developing so 
called Computational Hadronic Model (CHM) as an extension of FeynArts and FormCalc packages. 

The main purpose of FeynArts is to generate graphical and analytical representations of an unevalu- 
ated amplitude for scattering or decay processes. The results are produced by using Feynman rules as 
specified in the model files and do not include any further simplifications like index contraction or tensor 
decomposition and reduction. FormCalc, although is works as a shell for Form program, can manage the 
explicit analytical evaluation of the one loop integrals employing dimensional regularization and tensor 
decomposition techniques. Later, the renormalization can be performed by using either On-Shell, Minimal 
Subtraction (MS) or Constrained Differential Renormalization (CDR) schemes giving final results free of 
ultraviolet divergences. At the final stage, an automatically generated FORTRAN code can be used for 
numerical calculations of cross section or decay widths. 

Let us start with a description of the fields participating in CHM. It is only naturally to form three 
classes of fields: the octet of baryons B (16 states in total with octet of pentaquarks included), the decuplet 
of resonances T'^ (20 states in total with anti-decuplet of pentaquark resonances included), and the octet 
of mesons P. The external wave function for baryon (spin J = 1/2) is given by the usual Dirac spinor, 
and for the decuplet of resonances with spin = 3/2± we use the Rarita-Schwinger (RS) field described 
by 

1 1/2 

{k, m)=Y, ^^j^j • ^ ^) (11) 

i=-li=-l/2 

where CCfl^'t, are Clebsh-Gordan coefficients defined as 



V3(1+J)K1-J)!(U^)!Q-^)!' ^ ^ 
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The index m can take values {3/2, 1/2, —1/2, —3/2} and describes the projections of the = 3/2^ spin 
state. In our model file, the field enters as a product of (k) ■u{k) , specifying neither the polarization 
nor spin state. The explicit structure of the decuplet field (see Eq ITT]) is accounted for later, when the 
helicity matrix elements are computed. For the polarization vectors, we employ the following basis in the 
center of mass reference frame: 



,Esm{e) ,0,Ecos {9)) /m 



e>^^ (k) = (0, - cos (6) , ^i, - sin (6)) /V2. 



(13) 



Here, the angle 9 determines a direction of k vector. The external wave function for the meson field P is 
equal to one. The spin 3/2 propagator of the field T'^ is given using convention defined in [14|. 

Since the definition of coupling matrices T in the FeynArts is based on the use of the chirality projectors 

1±75 
2 



as 



) 



(0) (1) 

9r 9r 

(0) (1) 

9l 9l 



(14) 



the couplings are introduced in two model files responsible for the generic and particle levels, respectively. 
The first part of Eq.([14| describes the kinematic part of the coupling derived from the chiral Lagrangians 
given by Eqs.([Tl [7]and[9|). The second part of the Eq. (fT4l) represents the coupling strengths defined at 
tree level (first column) plus counter terms (second column). Coupling strengths g^^^ in Eq. lfHl) are 
represented by the Clebsh-Gordan (CG) coefficients of SU (3) group normalized by pion decay constant 
/tt for the decays or production channels. For the hadronic interactions, constants g^^^ are given by 
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The indexes R and L correspond to " + " and " — " in the couplings, respectively. The fields Tj^fc 
Pij are given by Eqs.ffSl [8l [TO]) . but only numerical coefficients in front of these fields are used in Eq.(fT5 
For the QED decays we have 
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Here Qt,b,-k is a charge of {T^, P} fields given in the units of electron's charge. SU (3) electromagnetic 
charge matrix Q is defined as 

2/3 

Q=\ -1/3 I . (17) 
-1/3 
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Figure 1: Next-to-the-leading hadronic decay with the adjacency of order five for the decay Pt^+- 



The counter terms of equation E q.(|14D are given through the wave function renormahzation of the baryon 
field. Here we adopt notation of [la] 



B) 



& - B) 



where indices {A,B} correspond to the external baryons of the process A ^ B. The wave function 
renormahzation constants |<J/^ij, (^/^/jI are computed through the truncated self energy graphs and 
given as an example for the baryons of the type A 



Re 



, (19) 



where L, R and S are left-handed, right-handed and scalar parts of the truncated self energy, respectively. 
Another counter term arising from the radiative decays of the decuplet into the octet state {T — > i? + 7} 
is described by the Lagrangian 



£TB7 = iQ_^B-f^'J5QT''F. 



A 



(20) 



The unknown coupling constant has been determined in [13|] from the measured branching ratios of 
A Nj and H**^ — > 2^7. The Clebsh-Gordan (CG) coefficients for this counter term are simply given by 



9^^{ {T - B7) = ±Q^T''"'eadeQtBl, 



(21) 



and have dimension of 1/M. 

Due to the presence of the Svr interaction in the chiral theory, a coupling of the adjacency of order five 
(see, for example, FigH]) is introduced in the hadronic decays of the baryons. 

Since, the couplings of the order five adjacency (see Fig.([T])) described by the same CG coefficients as 
for the tree level couplings of order three adjacency, it is straightforward to program the same coupling 
constants for every pion field appearing in the loop. This can be easily accomplished if the field 11 is 
introduced in the loop. This field 11 appears only in the loops with the Stt coupling, with a propagator 
described as 



An(g)=^' 



1 



1=1 



m. 



2 ' 



(22) 



where rrii is the mass of the meson in the SU (3) octet (see Eq.(l3|)). The kinematic part of the coupling 
for this type of processes is exactly equal to the couplings of the hadronic interactions described by chiral 
Lagrangians (Eqs.([Tl [7] and [11)). The same is true for the CG coefficients, but with constant replaced 
by fl in the first line of Eq. (fT5]l . 
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IV. APPLICATIONS AND TESTS 

The extension towards the hadronic sector in FeynArts and FormCalc which we propose in this work 
can be appHed to the virtually any hadronic process. We choose two processes relevant to the problems 
of the hadronic physics. These are pentaquark photo-production and Compton scatte ring in the studies 



of the nucleon polarizabilities. Since these problems were addressed earlier by [18|| and [20(|, this will serve 



as a test of the rather involved calculations using our Computational Hadronic Model (CHM). 



A. Pentaquark Photo-Production 



The question of existence of exotic baryons is of crucial importance to the physics of the strong in- 
teractions. The bound states of five quarks (pentaquark) can be described by the at least two models. 
Nucleon-Kaon molecule model, proposed by |16[], assumes that pentaquark is formed by a nucleon and a 
kaon bounded via the nuclear strong interaction. An alternative model proposed by [171] treat a pentaquark 
as a system of two diquarks and one anti-quark {q^q), bounded via the colour strong interaction. Since 
the range of the nuclear strong interaction is of order of ~ 1 fm, compared to the color strong interaction 
range, of ~ 0.1 /m, the decay widths of the pentaquark in these two models are substantially different 
(~ 100 MeV vs ~ IMeV). Obviously, the detection and subsequent measurements of the pentaquark 
decay width will most probably rule out one of these models. In the experimental searches of the ex- 
otic baryons, attention goes towards lightest pentaquark state 0^ = {uudss}, which has a mass around 
1540 MeV . This state is attractive not only because it has the highest production probability, but also 
due to the fact that symmetries of the chiral Lagrangian forbid the electromagnetic decay of this state to 
the proton and photon at the tree level and hence makes state more stable with decay width probably 
around 1 MeV. Unfortunately, despite of tremendous experimental efforts, production of this state proved 
to be extremely problematic. The experimental difficulties emerge from the weak signal, low statistics, 
and high background in all production channels. 

Prom the point of view of relativistic ChPT it is possible to introduce pentaquark baryons in to the 
computational hadronic model due to the same symmetry of the chiral Lagrangian and similarities of the 
10 and 10 irreducible representations of the SU{3) group. If one assumes that octet of the pentaquarks is 
represented by the spin 1/2 states, and the anti-decuplet of resonances by the spin = 3/2^, it easy to 
construct a Lagrangian for both. Using notation of [2]|, we can describe the octet of exotic baryons by 
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and the anti-decuplet of pentaquark resonances by 
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The interaction Lagrangian for all the possible permutations of the fields B^T^,0^,T^^fj_ and can be 
written in the following form: 

£(1008) ^ Tr 057^75 {^M, O5} + 2^0 Tr 057^^75 [^^, O5] + THpf^Yl^A.T^n^^ + 

Cpo [O^A^T^T^ + f^A^T^Or,) + Cps (S^^FsTg'^ + f^A^T^B) + (25) 
2PoB (Tr 057^75 {^M' ^} + h.c) + 2^ob (Tr ^7/^75 [A^, O5] + /i-c) • 

Relevant to the physics of the pentaquark production, coupling constant Cpb in Eq. (f25l) is given by 
CpB = The constant qtep can be extracted from the information on hadronic decay width of 

the specific pentaquark state. The rest of the coupling constants can either be extracted using the Adler- 



Weisberger sum rules (see [2l|]) or require additional experimental input. The Clebsh-Gordan coefficients 
for the pentaquark states were computed in the same way as in Eq. ffTS]) and Eq. (fT6]l but with coupling 
constants taken from Lagrangian in Eq. (f25]l . 

It is natural to raise a question: if the does exist, how can chiral effective theory direct experimental 
searches for this state? One of the possibilities would be to look at the kinematic dependencies (see [l8| ) 
using our Computational Hadronic Model. There are three possible channels for the production 6^: photo- 
, hadro- and lepto-production. The most recent experimental efforts were directed to the photo-production 
of the 0^. Photo-production is described by the two reactions, 'yp —>■ Q^Kq on the proton target and 
jn Q^K- on the neutron target. The decision to look at these two channels explicitly is motivated 
by controversy over experimental results (see Our goal here is to construct the production cross 

sections for these two reactions and to compare our results to In the Born approximation diagrams, 
responsible for the production of are given by the contact, s, t and u channels and shown on the 
Fig. ([2]). Obviously there is a distinctive asymmetry in the production channels for the two targets. This 
asymmetry is reflected by the presence of the contact term in photoproduction on the neutron target. 



This result was pointed out originally by [l8| and confirmed in our CHM. According to the Lagrangian 



given by Eq. (f25]) . we can construct two parity states of 0^ with = 3/2"'" or = 3/2 with values for 



the coupling constant qknb = 0.53 for 0^(3/2+) and qkno = 4.22 for 0^(3/2") 18|| if we assume that 
this state has a narrow decay width around 1 MeV. The production cross sections for the contact, s— , t— 
and u- channels are shown on Fig.jH]) as functions of the photon's energy in the laboratory reference 
frame. The largest contribution comes from the contact term which is only present for the neutron target. 
Comparing our results to [l8| we can see what s- and u-channel cross sections are somewhat larger in our 
case. This can be explained by the fact that, unlike in CHM model, the couplings in [3| are adjusted by the 
four- dimensional formfactors, and, since in these channels baryons are further off-shell, the contribution of 
these channels is strongly suppressed. As for the contact and t-channel contribution, our results are close 
to (l8|l- The evident asymmetry between neutron and proton targets favors deuterium target (with proton 
as a spectator) for the photo-production of the 05"(3/2^) with production cross section around ~ 200 n6 
and ~ 30 nb, respectively. 



8 




Figure 3: The total cross section for the 0g photo-production in the contact, s-, t- and u-channels as a function of 
energy of the photon in the laboratory reference frame. 



B. Compton Scattering 



Although the CHM predictions of the pentaquark photoproduction confirmed the earher theoretical 
findings, our calculations are done using the Born approximation only, and hence it is reasonable to raise 
the question about the validity of our CHM extension of FeynArts and FormCalc for the Next-to-the- 
Leading-Order (NLO) calculations. 

To test our NLO calculations (not including resonances) we decided to compute amplitudes for the 
Compton scattering off the proton and neutron. If we consider a nucleon with structure, the Compton 
scattering amplitude in the rest frame has the following form fgcj]: 

e^Z^ -> 

MhN -/N) = ~e' ■ ~e + auj'u;~e' ■ ~e + pCe' x k')re x k ) + Oiuo^). (26) 

47rm 

Here, ("?, oj, k ) is the polarization vector, frequency and momenta of the incoming photon, respectively. 
Primed quantities denote the outgoing photon. Two structure constants a and j3 are the electric and 
magnetic polarizabilities of the nucleon. For the point like nucleon, we have (see Fig.(l4|): 



This amplitude, in order to satisfy the gauge invariance of the electromagnetic current {d^J^ = 0), is not 
renormalized, and expressed through the set of physical observables such as charge and mass. Moreover, 
in the Thomson limit when k' ^ k , the Compton amplitude M in Eq. (l26p should take the form of Mq 
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Figure 4: FeynArts output for the Compton scattering at the tree level in CHM. 




Figure 5: Representative diagrams for the Next-to-the-Leading order Compton scattering in CHM. 

in Eq. (f27l) according to the well known low-energy theorem. Hence, when soft photons are considered, 
the nucleon behaves as a point-like particle and structure constants of Eg. ((26]) could be calculated only 
through the loop diagrams. Although generally the loop contribution is ultraviolet divergent and requires 
renormalization, this is not a case for the Compton scattering. If M{k' — > fc ) — > Mq, Mq is not 
renormalized and the limit /c ' — > k does not remove ultraviolet divergences, then the amplitude M 
should be finite. If we take only the SU{2) triplet of the mesons, we will have to evaluate 22 diagrams for 
the neutron and 52 diagrams for the proton (see Fig. ([5])). It is necessary to mention that the diagrams on 
Fig.(l5]) represent only generic types of the topologies allowed in Compton scattering. When one employs 
CHM, the full set of graphs will have to include the crossed diagrams and wave function renormalization 
graphs absorbed into counterterms. Also, in the CHM, charged mesons are treated as non-selfconjugate 
fields which effectively induce double counting, so part of the final amplitude with the charged mesons in 
the loops has to be divided by the factor of two. 

Individually, the diagrams on Fig.(l5|) give the divergent contributions but their sum should be a finite 
result. That was a case when Compton amplitude was computed in the CHM. We observed analytically 
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that the final Compton amplitude is, in fact, finite (i.e. free of ultraviolet divergences). Moreover, we 
have successfully tested finiteness of the Compton amplitude for the SU{3) octet of mesons. Here we have 
calculated 44 graphs for the neutron and 104 (not including wave function renormalization graphs) for the 
proton and observed no ultraviolet divergences in the final amplitude. In the FormCalc, the final results 
(amplitude or Compton tensor) were presented analytically using the Passarino-Veltman basis. Due to the 
cumbersome nature of this type of calculations we leave analytical details out of this article. Although we 
can proceed and express polarizabilities of the nucleon numerically, we leave that for the next publication. 
We have satisfied purpose of this article already when ultraviolet finite results for the Compton amplitude 
were obtained up to NLO using our computational hadronic model. 



V. CONCLUSION AND OUTLOOK 

In the present work we have developed an extension (the Computational Hadronic Model) of FeynArts 
and FormCalc to include the hadronic sector using Chiral Perturbation Theory. We have included octet 
of mesons, baryons (and pentaquark baryons) plus the decuplet of resonances (and the anti-decuplet of 
pentaquark resonances). The CHMs provides a robust prediction of the asymmetry in pentaquark 
photoproduction between neutron and proton targets. It can be concluded that the asymmetry is caused 
by the presence of the contact term in the 7n — > -ftr~05' channel. 

The ability of the CHM to handle Next-to-the-Leading-Order calculations was tested through the fact 
that the Compton scattering amplitude should be ultraviolet finite and should not be renormalized. As 
was expected, even when the entire octet of mesons was included, the final amplitude did not exhibit any 
ultraviolet divergent behavior. The NLO results of the CHM can be expressed analytically using FormCalc, 
with an amplitude expressed in the Passarino-Veltman basis as well as numerically with application of the 
package LoopTools. 

Future applications of the CHM can be envisioned for calculations of many processes, including the full 
kinematic dependencies for the production/decay channels of hadronic physics. Source files of the CHM 
are available by request from the authors. 
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